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ABSTRACT 


This report presents a discussion of mathematical models of Inso- 
lation characteristics suitable for use In analysis of solar energy 
systems and shows why such models are essential for solar energy system 
design. A model of solar radiation fcr the Mojave Desert area Is pre- 
sented with probabilistic and deterministic components which reflect 
the occurrence and density of clouds and haze, and mimic their effects 
on both direct and Indirect radiation. The model has the capability of 
producing any or all of the following outputfK 

(1) A "clear sky" theoretical amount of solar radiation. 

(2) Solar radiation for a clear sky or a cloudy sky or for a sky 
partially clear and partially cloudy depending on certain 
probabilistic parameters. 

(3) An array of average solar energy reception rates (solar In- 
tensities) In kilowatts per square meter (kW/m2) for a 
specified length of time. 

Multiple comparisons were made between measured total energy received 
per day and the corresponding simulated totals. The simulated totals 
were all within 111 of the measured total. The conclusion Is that a 
useful probabilistic model of solar radiation for the Goldstone, Cali- 
fornia, area of the Mojave Desert has been constructed. 
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SECTION I 
INTRODUCTION 


After many years of study and experimentation, and in the face of 
continuing uncertainty over supply and price of conventional energy 
sources, serious consideration is now being given to the question of 
using solar-powered energy systems on a relatively large scale. A 
necessary precursor to construction of well-designed, efficient, and 
economically viable solar energy systems is the engineering analysis not 
only of the systems themselves but also of the solar radiation tnat will 
drive them. This report presents a discussion of mathematical models of 
insolation characteristics suitable for use in analysis of solar energy 
systems and shows why such models are essential for solar energy system 
design. The currently most-used Insolation model is described, as well 
as the improvements that might be made in it to suit it better for use 
in designing solar energy systems. The design and construction of an 
upgraded model is presented together with its preliminary accuracy 
testing . 
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SECTION II 

WHERE AN INSOLATION MODEL IS USED 


Before construction of any energy system is undertaken, there must 
be reasonable assurance that it will meet the demand it was planned to 
satisfy, and that it will do so with a low enough life-cycle cost to 
make the project economically attractive. System performance models can 
be used first to Judge design alternatives against each other and 
against criteria for performance and cost and then to alter the design 
of the most promising systems to improve performance and/or lower costs. 
The function of an insolation model can be better understood by looking 
at some of the essential features of a solar energy system model. 

Solar-powered systems can cover a range of applications, from 
space and water heating in a single structure to central station 
generation of electricity, and can vary widely in complexity. A 
generalized solarthermal system will be made up of solar collectors, 
possibly some sort of storage subsystem, and a subsystem to convert 
thermal energy to the desired form. Depending on the application, each 
of the subsystems might be quite simple or very complex; for 
Illustrative purposes it will be sufficient to think in terms of the 
general groups. In a gross sense, the energy output of a solar- powered 
system is determined over a given time period by the amount of solar 
radiation collected by the system and the overall system efficiency. 
The efficiency with which the system operates depends in turn on 
characteristics of the included subsystems and the parameters on which 
their individual performances will depend. 

Of the solar energy hitting a collector, a fraction, depending on 
the sun's position relative to the collector surface and the collector's 
own geometry and optical properties, falls to get to the absorbing sur- 
face. A portion of the absorbed energy is lost via heat leaks and 
reradiatlon; the amount is determined by collector properties and the 
temperature at which it operates along with other factors like ambient 
temperature end perhaps wind speed. The remaining energy is transferred 
from the collector as sensible heat in a fluid at a temperature 
depending on fluid characteristics, the temperature of fluid entering 
the collector, and the collector temperature. The temperature of fluid 
entering the collector depends on how much heat is removed from it by 
other subsystems, such as that devoted to energy conversion. The amount 
of heat required by the conversion subsystem is governed by the load it 
is to satisfy, by its own Internal properties, by the temperature of the 
heat supplied to it, and by the temperature of the sink to which it 
rejects heat (if it must). The characteristics of the storage subsystem 
exert an effect on both collector and conversion components. All of 
these influences are reflected in a subsystem model that is made up of 
an Interrelated set of mathematical models representing the performance 
of each component. Because each piece depends strongly on factors that 
vary significantly with time, the resulting model should reflect the 
Important dynamics. Inaccuracies associated with the various component 
models will propagate and compound during analysis of the system, of 
course. This means that each of them must represent the performance of 


2-1 


78-11 


the associated subsystem with greater accuracy th.jn Is required of the 
whole system model. 

Given that a system model can be developed that will allow 
calculation of system output as a function of Insolation and other 
weather parameters, where are we then? A viable energy system, solar or 
otherwise, must be capable of supplying the output expected of it over 
the course of its useful lifetime. Conventional systems can be designed 
with the appropriate capacity and then provided with the amount of fuel 
necessary to do their job. Fuel for a solar-powered system, sunlight, 
is completely outside the control of man. Design of a solar energy 
system — the relative sizing and performance specifications of compo- 
nents — must be done then on the basis not only of its Intended output 
but also on the basis of the energy input that can be expected during 
its lifetime. The question arises — how does one supply appropriate 
values for the prime driving function, solar radiation, to allow an 
estimate to be made of system performance over the span of 10 to 15 
future years? That performance analysis must be accurate enough to per- 
mit design of a system that meets output criteria and cost criteria in a 
situation where compensation for even moderate uncertainty by oversizing 
components can be prohibitively costly. 

Experimental measurements of solar radiation intensity could be 
used to drive a system model. Such measurements are scarce, limited to 
a few locations, and more often than not of questionable accuracy 
(Ref. 1). Only rarely, in fact, have the needed aspects of incident 
radiation been measured. Elnpirical data suffer from a more fundamental 
deficiency, however. Using radiation measurements as input for a well- 
conceived system analysis may give a good estimate of how the system 
would have performed during the time the data were taken, but that 
estimate would only be good for the period in which the measurements 
were recorded. To arrive at the desired system performance it would be 
necessary, in addition, to simulate the system's behavior over that 
whole period, a procedure that could be unnecessarily costly and time- 
consuming . 

What is needed is a representation of Insolation characteristics 
that depicts those aspects of both its long-term and short-term behavior 
on which system performance depends, expressed in terms of a one year 
description. That one-year description may never match insolation 
behavior for a particular measured year, but would be extrapolatable to 
match closely all important aspects of insolation Integrated over a long 
time. The representation would, in short, be the output of a mathemati- 
cal model describing solar radiation. Such a model, along with a suit- 
able data base, would allow average or representative future behavior to 
be predicted, along with estimates of the frequency and magnitude of 
deviations from that average. As noted above, the accuracy of the out- 
puts from this model must be greater than the accuracy required of the 
outputs of the composite system model. 
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SECTION III 
THE ASHRAE MODEL 


An insolation model currently enjoying wide use Is that developed 
by the American Society of Heating, Refrigeration and Air-Condltloning 
Qigineers (ASHRAE) (Ref. 2). It was not created for the purpose of 
analyzing performance In solar-powered energy systems. Rather, Its 
Intended application was In estimating heat load on buildings for the 
purpose of specifying heating-cooling systems for installation there. 
The form of the ASHRAE model was dictated by its purpose — this Is the 
case with all mathematical models. In general, the effect of Insolation 
on cooling system requirements Is felt on clear sunny days; they specify 
the conditions with which a cooling system must cope. Only clear days 
are modeled by the ASHRAE equations. 

Some discussion of what happens to sunlight on Its way to Earth's 
surface will aid In dissecting the ASHRAE model. Energy emitted by a 
point on the sun arrives at the edge of our atmosphere In parallel rays. 
Its Intensity at that point depends on Earth's distance from the sun, 
and varies slightly with time of year. As the solar radiation passes 
through the atmosphere. Its direct normal Intensity (Intensity on the 
plane perpendicular to the ray bundle's direction) Is attenuated. Some 
of the energy Is absorbed by molecules of atmospheric constituents and 
some is figuratively knocked out of the bundle of parallel rays by 
molecular and particulate scattering. The degree of attenuation from 
these effects Is a function of the distance the radiation has to tra- 
verse In the atmosphere, and the concentration of absorbing and 
scattering species contained there. The ASHRAE model uses the following 
equation to mimic these Influences under clear day conditions: 

IdN = NAe-B/si'^ (1 ) 

I[)jj Is the direct normal Intensity of radiation at Earth's surface. N 
is a clearness ..umber that varies up or down slightly from a value of 1, 
depending on geographical location and season, and reflects the Inevit- 
able variation In clarity of what Is considered to be a clear day. The 
parameter A Is classified as apparent radiation at atmosphere's edge; It 
has a different value lor each month and Includes the combined Influence 
of the sun's distance from Earth and some atmospheric attenuation. The 
value of B, the atmospheric extinction coefficient, also varies monthly, 
reflecting the concentration of absorbing and scattering species. Sets 
of values for both A and B were determined by empirical curve fitting. 
That Is, they were the values that, when Inserted Into Eq. (1), produced 
values of Ip|^ that best matched data actually measured over a long 
period at a site with a defined clearness number of 1. Finally, 1/sln ^5 
(where Is the sun's elevation angle) approximates the distance that 
the parallel bundle of rays travels In the atmosphere, which varies with 
time of day and time of year. 

A unit area of surface at ground level will receive direct radia- 
tion I[) at a rate corresponding to the direct normal Intensity modified 
by the cosine of a, the angle between the direction of Incoming rays and 
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the direction perpendicular tu the surface. 


Id = IriN 

In addJt.ion to the energy arriving in a direct line from the sun, 
the surface in question will receive radiation from two other sources. 
Some of the scattered rays will, after bouncing about in the atmosphere, 
rsach ground level and the receiving surface, coming from all 
directions. Light th"t has been reflected from the surroundings will 
also be picked up. These two effects are treated in the terms 

^DS = CIdnFsS (3) 


and 


IdG = rloFsc 

where Ips I® Intensity of diffuse radiation coming from the sky, C is an 
empirically determined factor showing monthly variation, and Fss is a 
geometrical factor relating to the amount of sky in a position to radi- 
ate to the surface. IpQ represents radiation reflected onto the surface 
from the ground around it; Ig is the total radiation intensity falling 
on the ground (determined as for any surface); r is ground reflectance, 
and Fsg i® another geometrical factor. Analogous terms dealing with 
reflection from other surfaces might be required if the surroundings 
warranted. To sum up, the total radiation intensity received by a 
surface near ground level , according to the ASHRAE model , is expressed 
as 


I? = IdN cos a ♦ CIdnFcs ♦ rlGFgG (5) 

When the quantities included are properly evaluated, this model provides 
good approximations for total radiation intensity as a function of time 
during clear, that is, cloudless and hazeless, weather. This is one 
aspect of several needed for accurate analysis of solar system perfor- 
mance . 
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SECTION IV 

INTERACTION OF SOLAR RADIATION AND COLLECTORS 


There are many ways cf collecting solar energy. As far as their 
dependence on the characteristics of insolation is concerned, they may 
be classified in terms of the degree of concentration they involve. 
While solar energy is intrinsically of high quality, it arrives at 
Earth's surface widely distributed and must be reconcentrated to be put 
to useful work. Flat plate collectors use large areas of absorbing 
material to Intercept the radiation as it falls unaltered on the col- 
lecting surface. The resulting energy is removed as heat by a fluid 
circulating over the surface. A flat pla'^e collector can use all the 
radiation that hits it, but at high operating temperatures heat losses 
from the large area of hot surface limit its efficiency. Collector 
designs that concentrate the radiation before it strikes the absorbir.g 
surface seek to reduce these heat losses by cutting down on the surface 
area of hot material, allowing higher efficiency. In effect, the con- 
centrating collectors fbcus the light on a small absorbing area, from 
which heat is removed by a circulating fluid. Concentration ratio is a 
measure of the area over which radiation is captured relative to the 
area on which it is focused; the higher the concentration ratio the more 
precise focusing is required. 

Only direct radiation is useful to a concentrating collector. The 
ASHRAE model allows estimation of direct radiation on clear days, where 
it comprises a large fraction of the total incident light. Energy 
systems must also work on days that are not entirely clear. Then the 
proportion of diffuse radiation is much larger, and concentrating 
collectors will experience their own degradation of output. In 
comparing systems, one must decide whether large areas of possibly less 
expensive nonfocusing collectors, with high heat losses but capable of 
using all the components of Incident light, are more or less effective 
than perhaps smaller arrays of more expensive focusing collectors tha*^ 
will attain high temperatures more efficiently but can't use all the 
light. This comparison cannct be made without knowing the availability 
of both direct and diffuse radiation as a function of time for all kinds 
of days. 

The ASHRAE equations embody a semi-empirical, deterministic model. 
By dealing exclusively with one kind of day, a type that is practically 
eventless except for the rise and fall of the sun, they can be success- 
fully applied. A requirement for dealing with all days demands a model 
with probabilistic components as well as deterministic ones. An ideal 
model will reflect the occurrence and density of clouds and haze, and 
will mimic their effects on both direct and indirect radiation. For 
general application to all collectors, another phenomenon must be con- 
sidered. That is the circumsolar radiation. This is caused mainly by 
scattering of the sun's rays by Earth's atmosphere, and possibly also by 
refraction to a smaller extent, and is always present. On clear days, 
the effect is small and is limited to a narrow angular diameter about 
the sun's disk. On hazy days the turbidity of the atmosphere Increases 
the circ imsolar radiation at the expense of the direct component. It 
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also increases the angular extent of circumsolar radiation, which then 
merges with the diffuse radiation. The result is that concentrating 
collectors cannot focus the sun's disk sharply. Loss from atmospheric 
defocusing of the sun's image becomes more severe as concentration ratio 
increases. Neither this effect nor the frequency and nature of unciear- 
ness can be modeled deterministically at this time. Random variables 
must be employed to estimate their influence. 

Since the ASHRAE model performs well in predicting radiation on 
clear days, it forms a sensible starting point for first attempts at 
constructing a generalized model. We will concentrate on the terms 
described by Eqs. (1) and (3)f regarding the dear-day expressions for 
direct and diffuse radiation as being a baseline condition that is modi- 
fied by the random effect of the weather. Modifications would occur via 
the insertion of a pair of reindom variables (call them M and ra) , one in 
each equation. 


IpN = MNAe-B/sin f3 (6) 

Ids - “CIdn^ss 

Any terms for reflected radiation (Eo. 4) that might be required 
to model a situation will automatically be modified, since they would be 
derived from modified estimates of total radiation on the reflecting 
surface. On a clear day, both M and m would carry values of 1, and the 
original ASHRAE equations would stand. As "unclearness" increases, the 
value of M would vary on a short time scale — say, hourly. An additional 
variable could be inserted into Eq. (6) to model the circumsolar radia- 
tion, giving in the end 


Id = DMNAe“®/si" (8) 

An estimate of how the variable D mignt depend on concentration ratio 
has been made, but must be verified. Equations (6) and (8) form a solar 
radiation model that should be much more suitable than either experi- 
mental measurements or the ASHRAE model for supporting engineering 
analysis of solar energy systems. The improved model will be more re- 
presentative of solar behavior than a set of measurements, and will deal 
somehow with all the aspects of radiation that are important to a 
collector . 

The question arises — where do values of M, m, and D come from? 
Their oasis is a series of simultaneous measurements made as accurately 
as possible in one location over as long a period of time as practica- 
ble. The parameters collected would include total insolation intensity 
and direct radiation as measured by a set of devices with a number of 
different concentration ratios. Data would be taken at small time 
intervals, on the order of a few minutes. Diffuse radiation intensity 
could be derived by differencing measurements of total intensity and 
direct radiation determined without concentration. These data and 
Eqs. (6) and (8) would be used to calculate simultaneous values for the 
three random variables. From the calculated values, a Joint probability 
density function for M and m would be determined, as well as a 
functional relationship between D and concentration ratio. The derived 
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probability density functions are used, while model calculations are 
being carried out, to generate values of the random variables which will 
lead to a time series of calculated Intensities with the same 
statistical properties as the original measurements. 

For probability density functions to give the most representative 
results, they should be based on a very large number of measurements. 
It is often the case that their form can be determined from a more 
limited set of data, after which they can be upgraded by small adjust- 
ments in the parameters in which they are expressed as more and more 
data become available. Modification of the density functions for appli- 
cation to another location may bt' possible, using a limited number of 
insolation measurements at the new site and correlation with other 
weather data that might be more abundant. Verification of such a trans- 
fer, and establishment of the conditions under which it would be valid, 
would require careful measurements for comparison with predicted values. 

A program for gathering meteorological and solar data using abso- 
lute calibration standards has been underway at the Goldstone Splice 
Communications Complex since June 1d7H (Refs. 3i '•). These measure- 
ments provide an archive of solar data with calibration traceable to the 
absolute reference scale established by the Kendall Frlraary Absolute 
Cavity Radiometer (PACRAP) II at the World Meteorological Organization 
Radiometer InterciMnparison at Pavos, Switzerland, in 1Q7‘'. Irradlance 
measurements referred to this scale are found to be numerically 2.11 
greater than those same measurements referred to the WSb International 
Pyrhel iometric Scale (IPS-IQ^b). In addition, the Goldstone 
measurements provide a data base that has allowed the development of a 
preliminary solar model such as the one sketched, which is described in 
the following sect lens. It must be kept in mind, however, that the 
approach to an accurate insolation model described in these paragraphs 
is only a beginning. This probabilistic model is still very much 
simplified, and future effort might profitably bt* spent investigating 
those factors, now determined empirically, to more precisely identify 
and separate their deterministic and probabilistic components. 
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SECTION V 
THE SOLAR MODEL 


A. INTRODUCTION 


The mathematical model simulates solar radiation for the Goldstone 
area and has the capability of producing any or all of the following 
outputs : 

( 1 ) A '’clear sky" theoretical amouvit of direct and total solar 
radiation. 


(2) Direct solar radiation for a clear sky or a cloudy sky or 
for a sky partially clear and partially cloudy depending on 
certain probabilistic parameters described below. 

(3) An array of average hourly solar energy reception rates 
(solar Intensities) In kilowatts per square meter (kW/m2) 
for 1 and 2, above. 


The first two can be provided for any desired period of time from 
2 minutes to one year. Values are In kilowatt-hours per square meter 
(kWh/m2). Averages other than hourly can also be provided by a small 
modification to the program. If the time period selected Is more than a 
day then dally total Insolation values are also provided. 


The starting point Is the ASHRAE clear day model which has been 
described above. This clear day model Is modulated by the effects of 
clouds. The distribution of clouds for any given time Is determined by 
the combination of statistical procedures, measured Insolation values 
over a six-month period and a data bank of 19 years of cloud cover 
information. The method is described in the following sections. 


B. MODEL DESCRIPTION 

Figure 5-1 is a simplified overall flow-diagram of the model. Each 
block Is explained below. 


1. Inputs 

The Inputs are divided Into two kinds; user Inputs and built-in 
inputs. The user must supply information pertaining to the geographical 
location of interest; station altitude above sea level and station 
latitude. These data are required in the computation of theoretical 
clear day insolation levels. The user must also supply the time that he 
or she is interested in (a single point) or the time interval of 
interest. The built-in data consist of the theoretical clear day model 
(ASHRAE model), the measured pyrheliometer data and the weather data. 
The latter are 19 years of cloud cover and other meteorological data 
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Sinplified Flow Diagram of the Solar Model 
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recorded at Edwards Air Force Base, a few tens of kilometers from 
Goldstone. 


2. Sunrise and Sunset Times 

Sunrise and sunset times for the day(s) of interest are computed 
from the ASHRAE model . 


3. Expected Number of Dropouts 

If a cloud passes between the sun and a pyrheliometer aligned 
along the direct solar ray the detected insolation level may drop 
significantly from the clear day ’evel for that instant of time. When 
the solar intensity drops significemtly (5!l or more) from the clear day 
level this is termed a "dropout." Figures 5-2 and 5-3 are plots of the 
output from a pyrheliometer at Goldstone for two different days in 1975. 
Figure 5-2 shows a cloudless day and Figure 5-3 shows a day when clouds 
produced one dropout. 

The model computes an expected number of dropouts for the day of 
interest based on a histogram, made from the measured pyrheliometer 
data, of the number of dropouts that occurred per day. This number of 
dropouts is transformed into a continuous variable and is used as an 
input parameter for computing the expected time durations of dropouts 
and clear sky states. 


It is necessary that the expected number of dropouts be character- 
ized by a continuous function because it is used to generate the 
expectation parameter of the exponential distribution which is continu- 
ous. The discrete distribution of measured dropouts per day from the 
pyrheliometer data is shown in Figure 5-^. The conversion of this dis- 
crete distribution of dropouts per day to a continuous distribution is 
achieved in the following manner. First, a random selection of an 
Integral number of dropouts per day is made based on the probability 
distribution of Figure 5-4. Then, with the assumption that the proba- 
bility of experiencing i to (i-fl) dropouts in a single day is uniform, 
the integral number of dropouts, i, is converted to a random value in 
the interval i to (i-fl). To do this the model simply uses the random 
number previously generated which determined the interval from i to 
(i-fl) in question. The conversion therefore takes the form: 


XDOTS * 1 -f 


XUNI - Pi 
PU1 - Pi 


(9) 


where 


XDOTS = the expected number of dropouts per day 

i X integer part of the expected number of dropouts 
XUNI X a random number generated on (0,1) 

Pi X the cumulative probability of observing (ifl) dropouts per 
day 
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NuMH« or drokxms in a day 


Figure 5-4. Histogram of the Frequency Distribution of 
Dropouts Per Day From Measured Data 


4. Determination of Sky Condition at Sunrise 

The solar model steps sequentially through the day of interest 
from sun-ise to sunset and at each step a decision is made whether the 
direct component of insolation is clear or intercepted by a cloud, and 
if so how much attenuation is experienced. The time length of the step 
can be chosen by the user; 5-minute time increments are probably the 
moat convenient and are used in this report. This procedure is used 
even if the user Interested in only a single point in the day. In 
order to step sequentially through the day it is necessary to determine 
the sky condition at sunrise as an Initial condition. The weather data 
input information consists of 19 years of hourly cloud cover data 
reported as percentages. The percentage cloud cover at sunrise on the 
day of Interest is found from the weather data. The user has the option 
of choosing a specific year of the weather data, choosing one year of 
the 19 year set by a random process or using 19 (or other) years 
averaged data. In any case the hourly data are extrapolated to the 
exact moment of sunrise. Assume, for example, that decimal percent sky 
cover at the moment of sunrise on the day of Interest is C.x. The model 
then generates a random number from the uniform distribution (0,1). If 
this random number is less than or equal to 0.x then the day of interest 
starts in a dropout condition. If the random number is greater than 0.x 
then the day starts in the clear sky condition. 
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5. Simulated Distribution of Duration Times of Dropout and 

Clear Sky States 

The Flow Diagram of the Solar Model in Figure 5-1 shows that the 
next procedure is the generation of a simulated distribution of the time 
durations of dropout and clear sky states. Figure 5-5 is a detailed 
Flow Diagram of this procedure which is Iterated as shown in the figure 
for each time step throughout the day of interest. 

The first step in this procedure is the computation of expected 
time durations of dropout and clear sky states for a specified time of 
day, namely the start of each 5 minute time step. Figure 5-6 shows a 
dropout and a clear sky state symbolically as step functions as they 
might be recorded from the output from a pyrheliometer . An expected 
time length of a consecutive pair of dropout and clear sky states is 
calculated. This time interval is ETOTAL as shown in Figure 5-6 and is 
given by 

ETOTAL = (TDOWN - TUP)/XDOTS (10) 

Where TDOWN is sunset time, and TUP is sunrise time and XDOTS is the 
expected number of dropouts per day computed at the time step of inter- 
est. The time interval, ETOTAL, is partitioned into an expected dropout 
time Interval and an expected clear sky time interval , ETIMEDO emd 
ETIMES, respectively. The partitioning is based on the percentage sky 
cloud cover for the specified time of day from the weather data. Thus: 

ETIMEDO = ETOTAL • PCT (11) 

where PCT is the percentage cloud cover (in decimal form) at the speci- 
fied time. It follows that the time duration of the clear sky state, 
ETIMES, is given by ETOTAL - ETIMEDO. 

At this stage the model has computed expected times of duration of 
dropout emd clear sky states based on the expected number of dropouts 
(derived from measured data) and the percentage sky cover both of the 
specified time of day. These expected time lengths. ETIMEDO and ETIMES, 
are recomputed at each 5 minute time step throughout the sunrise to sun- 
set time interval. 

The next block in the Flow Diagram of Figure 5-5 shows that a 
decision is made whether to continue in the given state or whether to 
switch states. Both dropout and clear sky state distributions are 
characterized by the exponential density function. 

f(t) = ae-at (12) 

where t is time and Va represents the expectation value for the func- 
tion. The parameter 1/a takes on the value of ETIMEDO or ETIMES depend- 
ing on sky condition— dropout or clear sky, respectively. The exponen- 
tial density function was chosen because of its property of "lack of 
memory." That is, th« probability of a future event is Independent of 
what happened in the past. Thus the probability of a dropout continuing 
is independent of how long it has been going on, and similarly for 
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Figure 5-5. Flow Diagram of the Generation of a Simulated Distribution of 
Times of Duration of Dropout and Clear Sky States 
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Figure 5-6. Step Function Illustration of a Dropout 
Emd a Clear Sky State 


a clear sky condition. No findings in the analysis of the measured 
pyrhellometer data deemed any other distribution more appropriate for 
the characterization of the two states. Integrating the density function 
over time (the 5 minute time increment) yields the probability of an 
event occurring. If, for example, the sky condition is in a dropout 
state then the Integral represents the probability of switching to a 
clear sky state. Conversely, if the sky were clear it would represent 
the probability of changing to a dropout condition. These probabilities 
are obtained directly from the exponential distribution function 

F(t) = 1 - e-^'t t 1 0 (13) 

A value for PR, the probability of changing states, is found by substi- 
tuting into Eq. ( 13 ) values for ETIMEDO (if the sky is in a dropout con- 
dition) or ETIMES (if the sky is clear) and time, t, (5 minutes). The 
probability of not changing states is 1 - PR. Next, a random number 
from the unifonn distribution on (0,1), XUNI, is generated. If XUNI is 
less than or equal to 1 - PR then the sky remains in the given condi- 
tion, the time of day is advanced by one step of 5 minutes and the pro- 
cedure is iterated throughout the day. 

If XUNI is greater than 1 - PR then the decision is made to switch 
states. This takes place at some random time in the 5 minute interval. 
For convenience assume that the sky condition is in a di^pout state. 
The exact instant in the 5 minute interval that the sky changes from a 
dropout condition to clear is obtained from the Inverse function of the 
exponential distribution function, which is given by 

t = - [log(XUNI)] (1/a) (14) 


5-9 


78-11 


where 1 /q is the present value for ETIMEDO and t presents elapsed time 
in the 5 minute interval prior to switching states. In Eq. (14) t will 
always be given by 0 < t < 5 because a decision to change states cor- 
responds to a value of XUNI given by 1 - PR < XUNI < 1 . Therefore gen- 

erating a value for XUNI which leads to a decision to change states 
* results in the production of a random time within the 5 minute interval 
at which the switch is made. The value of t (in minutes) is subse- 
quently added to the time of day corresponding to the beginning of the 5 
minute interval auid the resulting time represents the end time of the 
dropout state and the beginning of the clear sky state. This procedure 
starting with the computation of a percent sky cover, PCT, and an 

expected time duration of the present state is repeated over the entire 
sunrise to sunset time, as shown in Figure 5-5. There are certain 
conditions which, if they exist, cause a deviation from this 

probabilistic method of deciding when the dropout and clear sky states 
begin and end. The conditions are: 

( 1 ) If the expected time length of whichever state we are in is 
greater than 0 and less than or equal to 0.001 then the 
exponent in the distribution function is sufficiently small 
that the function's value is undefined. Further, the 
probability of continuing in the state is essentially zero. 
Therefore, the state is automatically terminated at a random 
time within the five minute interval. 


(2) If the percent sky cover is 1.00 (lOOH) then the model 
automatically switches to or remains in the dropout state . 

(3) If the percent sky cover is 0 then the model automatically 
switches to or remains in the clear sky state. 


The last box of Figure 5-5 shows that the product of this 
procedure is an array of time durations, and times of occurrences, of 
dropout and clear sky states. 


6. Computation of an Array of Normalized Amplitudes of Dropouts 

The flow diagram of Figure 5-1 shows that the next step is the 
computation of an array of normalized amplitudes of the dropouts. The 
normalized amplitude of a dropout is the ratio of the measured energy, 
in this case the pyrheliometer data, to the clear sky theoretical energy 
derived from the ASHRAE model. Thus C ^ AMP ± 1, where AMP is the 
normalized amplitude of a dropout, and AMP = 1 for a clear sky. 

Figure 5-7 shows a flow diagram of the computation of an array of 
normalized amplitudes of dropouts. A linear regression equation, 
Eq. (15), of the amplitudes of dropouts on their lengths (time 
durat.'.ons) was computed from the measured pyrheliometer data. This 
regression equation is 

AMP s 0.44051 + 0.77653 (TD - 0.06155) (15) 
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where AMP is the normalized amplitude of the dropouts and TD is their 
time duration. Expected normalized amplitudes of the dropouts are com- 
puted from Eq. (15) using the distribution of dropout time durations 
found in Figure 5-5. 

There was a considerable amount of residual variation at each am- 
plitude level in the measured pyrheliometer data. This residual varia- 
tion was independent of the time duration of the dropouts and conse- 
quently was assumed constant. Furthermore, for a specified time dura- 
tion the amplitudes of the dropouts were assumed to be normally dis- 
tributed about the corresponding expected amplitude from the linear 
regression line. Therefore, after an expected amplitude level for each 
dropout is obtained from Eq (15) a randomized amplitude level for each 
dropout is generated from a normal distribution with mean AMP(i) and 
standard deviation SIGMA. AMP(i) represents the expected amplitude from 
Eq. (15) and SIGMA lf‘ the square root of the residual variation from the 
regression line. 

This process is Iterated, as shewn in Figure 5-7, to form an array 
of normalized dropout amplitudes for the entire sunrise to sunset time 
period or interval . 


7. Outputs 

At this pxDint the Solar Model has characterized the entire day 
from sunrise to sunset with a set of dropout and clear sky states, their 
time durations and times of occurrence, and a corresponding set of 
normalized amplitudes for the dropouts. This information along with the 
ASHRAE clear sky model can now be used to compute an array of quantities 
describing the amount of solar radiation available on any given day. 

It is known from experience with solar energy instrumentation and 
measurements at Goldstone that the clearness number in the ASHRAE clear 
day model, N in Eq. (1), varies from day to day and possibly even within 
a day. In order to Justify the use of the ASHRAE model, an evaluation 
of its accuracy with respect to the clearness number was made using the 
six months of measured pyrheliometer data from Goldstone. A random 
sample of thirty-five clear sky days (no dropouts) was taken from the 
measured data, and each of these day’s data was integrated over the 
entire day to yield a measured amount of energy received per day. In 
addition, for each day, a corresponding theoretical amount of energy was 
computed from the ASHRAE model using the value of 1.05 for the clearness 
number. For each pair of theoretical and measured quantities a ratio 
was formed as well as an average ratio and standard deviation. On the 
average the measured radiation was &. 9 % less than the theoretical radia- 
tion. This corresponds to an average measured to theoretically expected 
ratio of 0.9113312 with a standard deviation of 0.04964. The average 
value by which the theoretically expected radiation exceeded the mea- 
sured radiation can be related to the variability of the atmosphere at 
Goldstone. To compensate for this effect the Solar Model uses a clear- 
ness number N’, where N' comes from a normal distribution with mean 
0.95690 and standard deviation 0.04964, where 0.95690 = (1.05)0.9113312. 


5-12 


78-11 


A more accurate estimate of the appropriate clearness number of Gold- 
stone represents one area for future improvement of the model. 

For each successive time of day a clear sky theoretical solar 
intensity level (kW/m^), is computed from the ASHRAE model. This inten- 
sity is multiplied by a randomly generated clearness number N' from a 
normal distribution, as described above. Finally it is multiplied by 
the normalized amplitude value corresponding to the current time of day. 
This results in a solar intensity level for the specified time of day. 
The model treats each intensity level as if it were constant over a one 
minute time interval. By multiplying each ir tensity level by one minute 
the model computes the energy (kWh/m^) received each minute of the day. 
Integration over the day, or over any time interval (t^ - to) yields the 
solar energy received In that period. Figure 5-1 shows the various 
possible outputs. 
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SECTION VI 

EVALUATION OF THE MODEL AND CONCLUSIONS 


The frequency distribution of the number of dropouts per day from 
the measured data was presented in Figure 5-4. Figures 6-1 and 6-2 also 
present descriptive information on the measured data. Figure 6-1 is a 
frequency distribution of dropout times of duration computed from the 
measured data. Figure 6-2 shows the frequency distribution of the drop- 
out magnitudes, also computed from the measured data. A dropout magni- 
tude is defined as 

Magnitude = 1 - Normalized Amplitude. 

The first step in the evaluation of the Solar Model was to deter- 
mine if the model produced the appropriate number of dropouts as well as 
the correct distributions of amplitudes eind times of duration. The 
model was therefore evaluated against the measured data. Figures 6-3 
through 6-7 show the actual (from measured data) and four simulated dis- 
tributions of the dropout amplitudes as a function of their times of 
duration. It may be seen from the figures that there is a close simi- 
larity between the actual and simulated distributions. The total number 
of dropouts for the actual and the four simulated data sets were: 249, 
189, 200, 192, and 187 respectively. The only distinct difference 
between the actual and simulated distributions was that the actual dis- 
tribution exhibited a higher frequency of shorter dropouts than the 
simulated data. This explains the overall number of dropouts being 
fewer in the simulated data. 

A comparison was made between the actual and simulated frequency 
distributions of dropouts per day, as shown in Table 6-1. 


Table 6-1 . Time-'^f-Day Distribution of Dropouts for 
Actual Data and Four Simulations 


Before 10 a.m. 

10 a.m. -2 p.m. 

After 2 p.m. 

Actual 

94 

125 

143 

1st Simulation 

91 

117 

117 

2nd Simulation 

105 

124 

113 

3rd Simulation 

100 

104 

115 

4th Simulation 

97 

110 

118 
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Kluui'f* 6-1. Frequonoy Distribution of Dropout Timss of 
Duration From th<* Measured Data 


If any portion of a dropout was eontalnod in a time Interval, the 
dropout was olasslfled as occurring in that time perlc'd. This gives the 
dropouts the possibility of occurring in more than one time Interval and 
results in a larger exi>ected number of dropouts in the later time inter- 
vals of the day. The day was broken up into three time periods: before 
10 a.m., between 10 a.m. and 2 p.m., and after 2 p.m. The observed fre- 
quency of occurrence of dropouts during these three intervals was QU, 
1.?S, .and 143 respectively. The uniform Increase in the frequency of 
dropouts occurrlt\g in the later times of the day, as expected, coupled 
with the limited amount of measured vlata available, led to the decision 
not to weight the probabilities of dropouts occurring as a function of 
the time of day. This is another area wtiere future refinement of the 
mode] could prove uset\:l. 

Multiple oomixarisons were also made between measured total energy 
received per day and the corresponding simulated totals. The simulated 
totals vary because of the probabilistic nature of the model. Table 6-2 
is a summary of the actual and simulated data sets. The simulated 
totals were all within Ilf of the measured total. The solar model con- 
sistently provided total energies which were conservative estimates of 
the available energy. 
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Figure 6-2. Frequency Distribution of Dropout Magnitudes 
From the Measured Data 



Figure 6-3 • Bivariate Plot of Dropout Amplitude 
as a Function of Time of Duration 
for the Measured Data 


6-3 



NdWALIZCO AMPLITUDE _ NOPMALIZED AMPLITUDE 


78-11 



'igure 6-4. 


Bivariate Plot of Dropout Amplitude 
as a Function of Time of Duration 
for the First Simulation Run 



TIME OF DURATION (DECIMAL DAY) 

Figure 6-5. Bivariate Plot of Dropout Amplitude 
as a Function of Time of Duration 
for the Second Simulation Run 
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Figure 6-6. Bivariate Plot o'* Dropout Amplitude 
as a Function of Time of Duration 
for the Third Simulation Run 



IIME OF DURATION (DECIMAL DAY) 

Figure 6-7. Bivariate Plot of Dropout Amplitude 
as a Function of Time of Duration 
for the Fourth Simulation Run 
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Table 6-2. Actual and Simulated Total Energies. (Also Presented 
Are the Percentages These Totals Represent of the 
Actual . ) 


• 

Total Received 
Solar Energy 
kWh/m2 

% of 
Actual 

Actual 

733.37 

100 

1st Simulation 

675.53 

92.1 

2nd Simulation 

698.85 

95.3 

3rd Simulation 

661.92 

90.3 

4th Simulation 

658.78 

89.8 


Table 6-3. 

Percentages 
Correspond 
the ASHRAE 

to Which the Energy 
to the Expected Clear 
Model 

Totals (From Table 6-2) 
Sky Total Produced from 




% 

Clear Sky Total 

ASHRAE Clear 
Sky Total Energy 



100.00 


Actual 



72.1 

1st 

Simulation 



66.4 

2nd 

Simulation 



68.7 

3rd 

Simulation 



65.1 

4th 

Simulation 



64.8 


6-6 





78-11 


Finally, comparisons were made for the sampled days, between the 
ASHRAE clear sky expected solar energy, and the actual and simulated 
energy totals. Table 6-3 provides a sunmary of how these totals com- 
pared. The actual energy available observed over the 119 day period 
represented about 72H of the clear sky expected total, while, consistent 
with the findings cited in Table 6-2, the simulated energy totals 
represented slightly smaller percentages. The extreme divergence 
between the observed and expected clear sky total points to the need for 
a solar radiation model beyond the capability of the clear sky model for 
the design, assessment, or analysis of concentrating solar-thermal 
systems. 

Overall the model works well. It must be noted that the model is 
based on only 6 months of measured data. When more data are available, 
the model will be Improved even thougn it works well and is useful in 
its present state. The model produces estimates of solar energy which 
are conservative but this might be .unproved with further development. 

Several areas of refinement fbr the model were mentioned in the 
report . These are : 

(1) The use of more measured data on which to base the model and 
characterize the dropouts. 

(2) Incorporation of cloud layer data from the weather data in 
generating dropout amplitudes. 

(3) Investigation of the dependence of the probability of 
dropouts occurring as a function of time of day. 

(4) Refinement of the procedure for generating expected time 
lengths of both dropout and clear sky conditions, which is a 
critical procedure in the model. 

The conclusion is that a useful probabilistic model of solar 
radiation for the Goldstone area has been constructed. This is a pre- 
liminary model which works well and which might be Improved in several 
areas with further work. 
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